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Abstract. Let g be a scattering metric on a compact manifold X with boundary, i.e., a 
smooth metric giving the interior X° the structure of a complete Ricmannian manifold 
with asymptotically conic ends. An example is any compactly supported perturbation of 
the standard metric on M n . Consider the operator H = h A + V, where A is the positive 
Laplacian with respect to g and V is a smooth real-valued function on X vanishing to 
second order at dX . Assuming that g is non-trapping, we construct a global paramctrix 
U(z,w,t) for the kernel of the Schrodinger propagator U(t) = e~ ltH , where z,w £ X° 
and t 7^ 0. The parametrix is such that the difference between U and U is smooth and 
rapidly decreasing both as t — * and as z — + dX, uniformly for w on compact subsets of 
X° . Let r = a; -1 , where x is a boundary defining function for X, be an asymptotic radial 
variable, and let W(t) be the kernel e~ tr / 2t U(t). Using the parametrix, we show that 
W(t) belongs to a class of 'Legendre distributions' on X X X° xR>q previously considered 
by Hassell-Vasy. When the metric is trapping, then the paramctrix construction goes 
through microlocally in the non-trapping part of the phase space. 

We apply this result to obtain a microlocal characterization of the singularities of 
U(t)f, for any tempered distribution / and any fixed t j£ 0, in terms of the oscillation of 
/ near dX. If the metric is non-trapping, then we obtain a complete characterization; 
more generally we need to assume that / is microsupported in the nontrapping part of 
the phase space. This generalizes results of Craig-Kappeler-Strauss and Wunsch. 

1. Introduction 

Let (X, g) be a scattering manifold of dimension n. Thus, X is a compact n-dimensional 
manifold with boundary, and g is a metric in the interior of X taking the form 

dx 2 h 
■9 = — + — 

near the boundary. Here a; is a boundary defining function, y are local coordinates on the 
boundary Y = dX extended to X and h is a 2-cotensor that restricts to a metric (i.e., is 
positive definite) on Y . We shall assume that a; is a globally defined, smooth function on X 
that vanishes only at dX, and denote r = x^ 1 , which is analogous to the radial variable on 
Euclidean space. The metric takes a more familiar form when written in terms of r: 

g = dr 2 + r 2 h(r~ x ,y, dr, dy), 

where h is smooth in the first variable and y (smoothness in i = 1/r £ [0, e) is of course 
a much stronger condition than smoothness in r £ (e _1 ,oo)). In fact, one can choose 
coordinates locally so that g takes the form 

h „ t 2 i 2i / — i dx 2 h(x,y,dy) 
(1.1) g = dr +r h{r ,y,dy) = — r -\ = ; 
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here h is a metric in the y variables depending parametrically (and smoothly) on x (see |1(J] 
for a proof). Thus, g is asymptotically conic; it approaches the conic metric dr 2 + r 2 ho, 
where ho = h(0,y,dy), as r — ► oo. The boundary dX is then at geometric infinity, with 
each point of dX representing an asymptotic direction of geodesies. Euclidean space, with 
its standard metric or any compactly supported perturbation of it, is an example. We shall 
assume henceforth that coordinates have been chosen so that the representation for the 
metric holds. Let A denote the positive Laplacian with respect to g, and let H = + V, 
where V E x 2 C°°(X; R). Thus V is a short-range potential 1 . We consider the time-dependent 
Schrodinger equation 



where z E X, t E R, D t — — * an d / is a given distribution on X. Let U(t) — e~ ltH be 
the propagator for H. We wish to construct a parametrix U(t) for U(t) which captures all 
the singularities of U; in particular, we want U(t) — U(t) to have a kernel which vanishes 
rapidly both as t — ► and as we approach dX. However, to simplify the construction we 
shall let only one of the variables (z,w) E X 2 , say the left variable z, approach infinity, 
while we shall only require uniformity over compact subsets in the w variable. 

To state our main theorem, we choose a function (f> E C°°(X) which is zero in a neigh- 
bourhood of the boundary of X, we let r = x^ 1 as above, and we define the kernel 
W(t) = e- tr2 l 2t U{t). We remark that we regard U(t) and W(t) as acting on half-densities, 
so the kernels contain a Riemannian half-density factor {dgzdgwdt] 1 / 2 in each variable. Then 
our main result is 

Theorem 1.1. Assume that the metric g is non-trapping. Then the kernel <p(w)W(t) is a 
fibred- scattering Legendrian distribution on X x X° x R + of order (| , j) in the sense of 
associated to the Legendre submanifold defined in Lemma Vd '.5\ 

Remark. In the trapping case, one can still construct a parametrix microlocally in the non- 
trapping region; see Section |SJ 

In more prosaic language, Theorem 11.11 savs that the propagator is given by oscillatory 
integrals of certain rigidly prescribed forms. Near t = 0, and with both variables w, z away 
from the boundary of X, these take the form of a Legendrian distribution (see Melrose- 
Zworski [H) 



where the integral is over a compact set K C M. k . Here, $ and a are smooth in all their 
variables. Near the diagonal, $ is given by dist(w, z) 2 /2, as in the flat Euclidean case, no v 
variables are required, and there is no integral. Associated to $ is a Legendrian submanifold 
of T*X x T*X x R, here given explicitly by 



(this is a somewhat non-invariant description, for precise details see Sj3J). The function $ 
becomes non-smooth outside the injectivity radius; geometrically this corresponds to the 
Legendrian L becoming non-projectable (i.e. the projection from the Legendrian to the base 
X x X is no longer a diffeomorphism) . The Legendrian, however, remains perfectly smooth; 
it is no longer parametrized by a function just of w and z, but one needs extra variables v 

lr The case of gravitational long-range potentials in the sense of |7| can also be treated. 



(1.2) 



(D t + H)u{z, t) = 0, u(z, 0) = f{z) E C-°°(X), 




L = {{w, t z, C, t) I (z, C) = exp g (u>, £), t = — } 
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(precisely, one needs at least k extra variables locally if the kernel of the projection from 
the Legendrian to the base has dimension k). 

When z approaches the boundary of X, we use local coordinates (a;, y) as described above, 
and sometimes use r = x~ l , the asymptotic 'radial' variable. In this region the propagator 
takes a more complicated form. We show that W(t) is a finite sum of oscillatory integrals 
of the form 

(1.3) i-"/ 2 - fe / 2 (27r)- fe f e l ^ x - v - w - v ^ xt a(x,y,w,v,t)dv, 

J K 

where again ip and a are smooth in all their variables, and one integrates over a compact 
set K cM. k . Notice that ip depends on all variables apart from t. 

Next, we address the question of determining the wavefront set of u(-, t) at a fixed nonzero 
(say, positive) time t in terms of the wavefront set of the initial data /. We are particularly 
interested in the case of 'interior singularities', lying above a point in the interior of X. It is 
well known that equation (|1.2|) has infinite speed of propagation, so to answer this question 
we must look for singularities of / at the boundary of X. Dually, since U(t) = U(—t)*, 
we can consider interior singularities of / and determine the singularities of u(t) that they 
produce. It is known from |21| that if there is an interior singularity (w, rf) for /, then for 
all positive times t there are quadratic oscillations of frequency l/2t in u(-,t) at infinity 
in the asymptotic direction of the geodesic emanating from (w,rj). Roughly speaking, the 
quadratic oscillations look like e lr / 2t in a conic neighbourhood of the asymptotic direction 
y; more precisely, this result is phrased in terms of the 'quadratic scattering wavefront set' 
(see jZU). 

The limitation of the result of is that different singularities along a single geodesic, or 
along a different geodesic with the same asymptotic direction, produce identical quadratic 
oscillation, hence consideration of the quadratic wavefront set alone will not result in precise 
propagation results. To analyse the finer structure of u(-, t), we divide «,(•, t) by the explicit 
quadratic oscillatory factor e lr / 2t , and find that the resulting function has oscillations which 
are linear in r that contain the desired information on the location of the interior singularity. 
The presence of linear oscillations of a function / is measured by the scattering wavefront set 
WF SC (/) (see [E]), a closed subset of sc Tg X X. (This bundle, defined in jgj is the restriction 
to dX of an rt-dimensional bundle over X which is a compressed and scaled version of the 
cotangent bundle T*X). We show that the asymptotics of geodesic flow on X determine 
two contact transformations Sf and Sb, which we call the forward and backward sojourn 
relations, from the sets F, resp. B C S*X° consisting of points in the cosphere bundle 
which are not forward, resp. backward trapped under geodesic flow, to sc Tg X X. They are 
related by F = — B and S/(C) = — S&(— ()■ The definition of Sf in local coordinates is as 
follows: Let 7(3) be the arclength-parametrized geodesic emanating from (w,r/) £ S*X°, 
let y = (yi, . . . , y n -i) be local coordinates on dX and let (y, v, fi = (/ii, . . . , /i„-i)) be the 
induced coordinates on sc Tq X X as in (|2.1I) . Then Sf(w,n) = (yo,u, fi) iff 

yo = lim 7(s) e dX 

s — >+oo 

is the asymptotic direction of the geodesic, 



(1.4) 



v = lim (s — rOy(s))), r = — 

s— >+oo v X 
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is the 'sojourn time' of the geodesic, and 

(1.5) /ii = V/i y Urn (y - y(l(s)))j)/x 

< 4 s — >+oo 
3 

measures its angle of approach to dX. (That the point j/o £ dX should be thought of as the 
asymptotic direction of the geodesic can be seen easily in the special case when X is the radial 
compactification of R ra .) Note that the sojourn time, thought of as depending on a geodesic 
and a point along it, is closely related to the classical Busemann function of differential 
geometry; see for example |17j . chapter 1, section 2. We use the term 'sojourn relation' for 
Sf since the coordinate v is analogous to the sojourn time considered by Guillemin [Hj. 

Theorem 1.2. Suppose f £ C-°°(X), and let u(; t) = e~ UH f . Let C = (w,r)) G B C S*X° 
be an non-backward-trapped point of S*X° . Then for any fixed t > 0, 

(1.6) C € WF(u(-,i)) if and only if ±S b (Q £ WF sc (e ir2 /"/). 

Here the factor t^ 1 acts by scaling the fibre variables. Similarly, if £ £ F, then for t < 0, 

(1.7) C G WF(u(-,t)) if and only if ±S f (C) £ WF sc (e' r2 / 2i /). 

Remark. The condition on the right hand side of l|1.6fl or i|1.7|) is not manifestly coordinate 
invariant. The sojourn relation changes under a change of coordinates, but so does the 
scattering wavefront set of e lr ^ 2t f, in such a way that the condition as a whole is coordinate 
invariant. A manifestly invariant description may be given in terms of the affine bundle of 
Lemma 13.31 

Example. Consider the free Hamiltonian H = |A on W 1 . The propagator in this case is 

U(t) = {2nt)- n/2 e l \ z - w \ 2/2 \ whence W(t) = ( 27 rf)-"/2 e -^-Wt e *H 2 /2t_ 

It is not hard to check that W(t) satisfies the conclusion of Theorem 11.11 in this case. 
Consider initial data 

/ = ( 27 rT)-"/ 2 e- l l z -"'l 2 / 2T , T > 0. 

Then the solution at time t = T is a delta function centred at w, hence its wavefront set is 

{(w,z) | z G S 71 - 1 }. 

In the free case, the backward sojourn relation is given by 

Sb(w, i) = (y — —z, v — —w ■ z, fj, = w — (w ■ z)z) , 

and we easily check that (|1.6(l holds in this case. 

Most previous work on Schrodingcr parametrices has focused on the case of flat space 
with a potential perturbation, where the geometric situation is substantially simpler. Very 
detailed parametrix constructions in this setting have been made by Fujiwara 0], Zelditch 
|23j . Treves [201 anc ^ Yajima [221 • m the case of curved space, very little was known. 
Kapitanski-Safarov |12j have shown that on W 1 with a compactly-supported, nontrapping 
potential perturbation, the fundamental solution is smooth for t > 0, and have exhibited a 
parametrix modulo C°°(IR n ) Such a parametrix, however, is not sufficiently specified 

at infinity to yield results about smoothness of the solution of the general Cauchy problem 
at t > 0. 
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Regularity results for the Schrodinger propagator for non-flat metrics in the form of 
Strichartz estimates have been obtained recently by Staffilani-Tataru |18| and Burq-Gerard- 
Tzvetkov £Q. Staffilani and Tataru ^B] proved Strichartz estimates for e~ ltH f where H 
is the Laplacian of a C 2 , compactly supported, nontrapping perturbation of the standard 
metric on R™, using the FBI transform and Littlewood-Paley decompositions to handle the 
rough metric. Burq, Gerard and Tzvetkov obtained Strichartz estimates, with a loss of 
derivatives compared to the flat Euclidean case, for compact manifolds or perturbations of 
the Laplacian on M n , without any nontrapping assumption. 

Various authors have considered the question of determining the singularities of u = 
e ltH f m terms of /. The first results about microfocal smoothness of u(t) for t > 
and general initial data were those of Craig-Kappeler-Strauss [3], who showed that, on an 
asymptotically Euclidean space, decay of the initial data in a microfocal incoming cone 
yields microfocal smoothness along the whole pencil of geodesies emanating from that cone 
for all t > 0. This result was refined in [23 > where the second author showed that absence 
of the quadratic-scattering wavefront set (see Section allowed one to conclude absence 
at varying times and along varying pencils of geodesies. This approach, while it specified 
in terms of the Cauchy data when and in what direction singularities might appear in A°, 
failed to say anything about where they might land. These propagation results have been 
extended to the analytic category by Robbiano-Zuily |16j . 

Our Theorem 11.21 gives a complete solution to the propagation problem in the case when 
the metric g is nontrapping, and in general, a complete characterization of the singularities 
of u(t) in the non-backward-trapped set F for t > 0, and in the non- forward-trapped set 
B for t < 0. Our results imply those of [22 > ancl hence those of [3j 2 , since the hypothesis 
on / required in |5J for microfocal smoothness of u(t) along all geodesies emanating from 
y 6 dX implies that in fact (y, v, fi) WF sc e ,r / 2t / for all (v, fx). Hence Theorem II . 21 yields 
the main boundary to interior propagation result of |21| as a special case. 

We thank Richard Melrose, Andras Vasy and Steve Zelditch for useful conversations, 
and the Erwin Schrodinger Institute, the Mathematics Department at SUNY Stonybrook 
and the Mathematical Sciences Institute at the Australian National University for their 
hospitality. 



We recall the definition of various structures associated to manifolds with boundary and 
corners needed in this paper. For further details, see Melrose ^3, Melrose- Zworski |14| . 
Hassell-Vasy 00. 

2.1. Scattering structure and Legendrian distributions. Let A be a d-dimensional 
manifold with boundary, and let a; be a boundary defining function for A. We identify a 
collar neighborhood of dX with dX x [0, e) so that specifying j/i, . . . yd-i local coordinates 
in d X gives local coordinates (x, y) for X . The Lie algebra of scattering vector fields V SC (A) 
consists of vector fields of the form ax 2 d x + xbid Vi with a, b{ £ C°°(A). Such vector fields 
can be described as the set of C°° sections of a vector bundle SC TX. The dual of this bundle 
we denote SC T*A; sections of it are locally spanned over C°°(A) by dx/x 2 and dyi/x. Hence 
any point q £ SC T* A has a unique expression 



2. Contact structures and Legendrian distributions 



(2.1) 




although we require more decay of our potential than is assumed in 
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which yields local coordinates (x, y, v, [£) for SC T*X, v and // being linear on each fibre. We 
say that a half-density a on X is a scattering half-density if it is a smooth and non- vanishing 
section of the bundle (A n ( ac T*X)) 1 / 2 ; such sections have the form a\dxdy\ . . .dyd-i/x d+l \ 1 / 2 , 
where a G C°°{X) is nonvanishing. The restriction of SC T*X to the boundary of X is denoted 

We define C°°(X) = f] l ^, Q x l C 00 (X), with its natural Frechet topology and denote by 
C~°° (X) its topological dual. We sometimes refer to these as the space of Schwartz functions 
and the space of tempered distributions on X, by analogy with Euclidean space. 

Now we recall some facts about the scattering calculus sc q> m - l (X; SC £1 1//2 (A)) acting on 
half-densities, which is indexed by two orders (m,l), the interior order m (which for a 
differential operator is the order of the highest derivative that occurs) and the boundary 
order I. The half-density factor will be understood from now on, and dropped from notation. 
The space fy m > l (X) is the same as x l ^ mfi (X), and * m '°(X) is a 'microlocalization' of the 
scattering differential operators of order m on X; it contains in particular all mth order 
differential operators generated over C°°(X) by V sc (^0- Operators P G sc if> m > l (X) arc 
determined up to sc $ m,,+1 (X) by the boundary symbol p, which is a smooth function 
on the boundary of the scattering cotangent bundle, 8C Tg X X, and up to sc ^ m ^+ 2 (X) by 
the boundary symbol p together with the boundary subprincipal symbol, which is again a 
smooth function p su ^ on sc Tq X X . The operator P is said to be elliptic at q G sc Tg X X if 
p(q) 7^ 0. The scattering wavefront set WF sc (u) of a distributional half-density u is defined 
by the condition that q G SC T£ X X is not in WF sc (u) iff there is an A G SC ^°'°(X) such that 
A is elliptic at q and Au G C°°(X), the space of Schwartz functions on X°. The scattering 
wavefront set of u is always a closed subset of sc Tg X X. We also have a scale of Sobolev 
spaces H%(X), defined by u G H£(X) iff Vi • • • V m u G L 2 (X) for all V u ..., V m e V SC (X), 
or equivalently, if Pu G L 2 (X) for all P e ^ rn '°(X). (Here L 2 (X) is defined with respect to 
the Riemannian measure dg.) 

For the purposes of this paper, we will often take the manifold with boundary to be 
X° x X° x K^ (this space is not compact, but that is irrelevant here). In that case the 
boundary defining function is t and local coordinates on the boundary will be denoted (z, w), 
where z G K" is a local coordinate for the first factor and w G M™ is a local coordinate for 
the second. In this case we use coordinates (t, z, w, r, £, rj) where we write points q' G 
sc r(i° x X° x R^ ) 

(2-2) l' = rd(l) + ± Cl d -^ + ±v/-f. 

i=l j=l 

Returning to the general situation, there is a contact structure defined at the boundary 
sc Tq X X of SC T*X. It is defined by the contact one-form 

(2.3) X = "(x 2 d x ,-)\{x = 0}, 

where u> is the symplectic form on T*X (which is canonically isomorphic to SC T*X over X°). 
In local coordinates \ = S Vidyi — dv, so \ 1S clearly nondegenerate. A change of boundary 
defining function x — > ax changes x according to x — » a\, so the contact structure defined 
by x is completely natural. A Legendrian submanifold of sc Tg X X is defined, as usual, to be 
a smooth submanifold of maximal dimension, namely dim X — 1, such that the contact form 
X vanishes on it. Any Legendrian submanifold L has a local nondegenerate parametrization 
in a neighbourhood of any q G L. This, by definition, is a function ip(y,v), with v G K fc for 
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some k > 0, such that 

dip 

the differentials dy— — ), i = 1, . .., k are linearly independent whenever d v ij) = 0, 

(2.4) dVl 

lb 

and L — {(y, d( — )) \ d v ih = 0} locally near q. 

The simplest situation is when L is projectable in the sense that the projection (y, v, /i) — > y 
from SC T^ X X to dX restricts to a diffcomorphism from L to <9X. Then y is a coordinate 
on L, so f is given by a function ip(y) on L. In this case, no extra variables v are required, 
and L is given by {(y, d(ip/x))} locally. If the kernel of the differential of the projection 
(?/, is, fi) — > y restricted to L has dimension k then at least k extra variables are required to 
locally parametrize L near q. 

A (half-density) Legendrian distribution of order m associated to to L is a half-density 
ua, where a is a scattering half-density and u is a finite sum of terms u « + w 0i where 
u 6 C°°(X) and is given by an certain type of oscillatory integral associated to a local 
paramctrization of L: 

(2.5) u t = (2n)- k x m - k / 2+d / 4 [ e l ^ v - v)/x a(y,v 7 x)dv. 

Here a is a smooth function of x,y,v with compact support and ip is a nondegenerate 
paramctrization of i on the support of a. The set of half-density Legendrian distributions 
of order m associated with L is denoted I m (X, L; sc f2 1 ^ 2 (X)), or just I m (L) when the space 
X is understood (in this paper, Legendrian distributions will always be half-densities). 

Legendrian distributions have a well-defined symbol map a m taking values in smooth 
sections of a line bundle S^(L) over L. This bundle is given by S^(L) = sc n 1/2 (L) ® 
\N*(dX)\ m ~ d/4 ®M®E, where Vt 1 ' 2 ^) denotes the half-density bundle over L, N*X is 
the conormal bundle, M is the Maslov bundle and E is the bundle described in jjj. To 
define the symbol, we choose d— 1 functions Xj in (y, u)-space which together with d Vi th give 
local coordinates in (y, w)-space. Then {A^} are local coordinates on L via the identification 
(|2.4|) . The symbol is given, using the identification l|2.4|l and up to Maslov factors, by 

-1/2 

|dA| 1/2 \{dib/dv = 0}. 

The symbol of u £ I m (L) determines u modulo I m+1 (L). Consideration of how the symbol 
changes under changes of parametrization ?/>, changes of coordinates (x,y), and changes of 
coordinates A, show that the symbol lives in the bundle S^ m \L) above If L is locally 
projectable, then the situation simplifies. We may take coordinates A on L to be y, and the 
v variables are absent, the determinant factor above is 1 and the symbol becomes 

a(0,y)\dy\ 1/2 . 

The important bundle in the factorization of S^ m \L) here is | N*(dX)\ m ~ d ^ 4 ; in particular 
this tells us that if we change boundary defining function from x to xa~ l {y) then a symbol 
of order m changes by a factor a m ~ d l i . 

We now recall the symbol calculus for a scattering pseudodifferential operator P e 
SC \|/*>0(X) acting on a Legendrian distribution u <E I m {L). In fact, we only need to consider 
the case when the symbol p of P vanishes identically on L, and in view of our application 
in 21 we use coordinates (|2.2() appropriate to the manifold M — X° x X° x K^o- Then 
the Hamilton vector field of the function p (extended into the interior of SC T*M arbitrarily), 



(2.6) 



» = a(0,y,v) 



det 



d(\,d v ip) 
d(y,v) 
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Hp, vanishes to first order at the boundary of SC T*M, so we define the rescaled Hamilton 
vector field sc H p to be x~ x H p restricted to sc T* dM M. Then Pu e I m+1 (L) with symbol 

(2.7) a m+1 {Pu) = ( - iCs. Hp - i(l +m - ^ + p sub )a"» ® |<ft|. 

Here d = 2n + 1 is the dimension of M, and £ir p denotes the Lie derivative acting on half- 
densities. The operation of tensoring with \dt\ is simply the natural isomorphism between 
S [m] {L) and S [m+1] (L). 

2.2. Quadratic scattering structure. We only touch on this very briefly since we only 
make fleeting use of it in this paper. The quadratic-scattering Lie algebra is defined by 
V qsc = xV sc (X). It is the space of smooth sections of a vector bundle i sc TX which is locally 
spanned (near the boundary) by x 3 d x and x 2 d Vi . The dual bundle is denoted ^ SC T*X; 
sections of this bundle are locally spanned over C°°(X) by dx/x 3 and dyi/x 2 . The vector 
fields in V qS cPO he in a calculus of quadratic scattering pseudodifferential operators much 
like the scattering calculus, and with a boundary symbol and associated wavefront set lying 
in qsc Tg X X. It is in terms of this wavefront set that the propagation results of [21] are 
couched, but it will not directly concern us here. 

2.3. Manifolds with corners with fibred boundaries. Above, we have looked at X° x 
X° x M^o- However, in our analysis of the Schrodingcr propagator U(z,w,t), we wish to 
let the first variable z approach the boundary of X. Hence we need to study the manifold 

(2.8) X 2 = X x X° x R >0 . 

We denote the two boundary hypersurfaces of X 2 at t = and at x = by mf and sf, the 
'main face' and the 'side face', respectively, and we denote the corner mf n sf by K. Thus 
mf = X X X°, sf = dX X X° X R^ and K = dX x X° . We shall see that the geometry 
of this manifold arising from the structure of the Schrodinger operator D t + H is that of a 
manifold with corners with fibred boundaries; see and |S] for a discussion of the general 
situation. Here, we restrict ourselves to describing the particular situation of X 2 , which is 
considerably simpler. 

We first describe the relevant fibration. Define (f> : sf — > dX x X° by projection off the 
last factor: 

(2.9) <f>(y,w,t) = (y,w)edXxX°. 

We then define a Lie algebra of smooth vector fields on X? as follows: if V is a smooth 
vector field on X 2 , then 

(2.10) V E V if and only if V(xt) = 0(x 2 t 2 ), V is tangent to 

the fibration <j> at x = and V vanishes at t = 0. 

This Lie algebra is independent of the choice of boundary defining function x. Corresponding 
to this Lie algebra is a vector bundle ^TX 2 such that is canonically identified with the 
space of smooth sections of ^TX 2 . In local coordinates, this is spanned (near the boundary 
x = 0) by 

(2.11) xtd Vi , xtd Wi , t(tdt—xd x ), tx d x . 

Note that t 2 dt and txd x are not separately in ^V, only their difference is. The dual bundle 
^T*X 2 is spanned over C°°(X t 2 ) by the one-forms 

j i 1 ^ xdt + tdx dy% dwj, 
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Thus every point in ^T*X 2 near {x = 0} may be expressed 

(Jl and £ are used since we are reserving the symbols fi and £ for later use), which yields 
local coordinates (x,y,w,t,cr,~p,£,T) on ^T*X 2 . A fibred-scattering half-density is defined 
to be a section of f\ n+ (^T*X 2 ), n = dim A, which is smooth and nonvanishing. Near the 
boundary, such a half-density has the form 



dx dy dw dt 



j-2n+2 x 2n+l 



1/2 

, where a 6 C°°(X 2 ) is everywhere nonzero. 



On ^T^Xf, (that is, at t = 0), uj{xt 2 d tl -) is a contact form (i.e. nondegenerate) for 
x > 0, but it degenerates at x = 0. This is evident from its local coordinate expression 

u(xt 2 dt, •) C {t = 0} = da + xdr + JLdy + £dw. 

As described in jgj, although this is degenerate at ^T^X 2 , that is, at ^T^Xf n {x = 0}, 
there is a natural fibration /3 from ^T^Xf to a bundle 7 over if, with a nondegenerate 
contact form on 7 such that ui(xt 2 dt, ■) restricted to {t = 0, x = 0} is the lift of the contact 
form from 7 to ^T^X?. Here the bundle 7 is naturally isomorphic to sc Tg XxXo (X x A ), 
the contact form is the natural contact form (|2.3|l on sc Tg XxXO (X x A°) and the fibration 
/3 is given by 

(2.12) (3(Q,y,w,0,<r,~p,l,T) = (y,w,a,JI,l) 

in the coordinates above. 

A fibred-scattering Legendrian submanifold J of ^T*Xf, as defined in [5], is a submani- 
fold of ^T^X 2 which is Legendrian with respect to the contact form aj(xt 2 d t , ■) for x > 0, 
which meets the boundary of ^T^Xf at x = transversally, and such that the fibra- 
tion /?, restricted to dJ, is a diffeomorphism from dJ to a Legendrian submanifold G of 
sc Tg XxXO (A x A ) (with respect to the natural contact structure on sc Tg XxXO (X x A )). 

A local nondegenerate parametrization of J near q £ d J, is a function ijj(x, y, w, v), with 
v £ M fc for some fe ^ 0, such that, locally near q, 

dtp 

the differentials d(— — ), 1 ^5 i ^ fe are linearly independent 

(2.13) dVl 

and J = {{y,dA) \ d v ^ = 0}. 
xt 

It then follows that the function ^(0, y, w, v) parametrizes the Legendrian G. 

A Legendrian distributional half-density of order (m, r) associated to J is a half-density 
ua, where 5 is a fibred-scattering half-density and u is a finite sum of terms Y] , Uj + 
Uo + Uf, where uq is a Legendrian distribution of order m associated to the interior of J 
(where the structure is locally of scattering type, so this has already been defined), Uf is 
a fibred Legendrian distribution of order r associated to G, which is supported away from 
mf (defined below) , and Uj is given by an certain type of oscillatory integral associated to a 
local parametrization of J. Namely, Uj is given by 

(2.14) ( 2 ^)-fe i ™-fc J /2+(2n+l)/4 a ,r-fe/2+(2n-l)/4 f e **i(.*,V,*',*)/'t aj ( X) y )W)V)t ) fa 

J E 

Here v € K fci , E C is compact, dj is a smooth function of x,y,w,v with compact 
support and ip is a nondegenerate parametrization of L on the support of a. Similarly, Uf 
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is given by a sum of terms of the form 

(2.15) ( 27r )-V- fe /2+(2n-i)/4 [ e^^- v ' w ' v)/xt b.j{x,y,w,v,t)dv. 

Je 

with b smooth and supported away from t = 0. It is straightforward to check that (|2.14l) 
and H2.15fl are compatible. Note that, for fixed t > 0, the expression l|2.15[) is a Legendrian 
distribution of order r — 1/4 associated to the Legendre submanifold t -1 G, where t^ 1 acts 
by scaling in the fibre variables. 

Fibred-scattering Legendrian distributions of order (to, r) on L have a well-defined symbol 
map a m taking values in x r ~ over L, which is defined by continuity from the 

symbol map in the interior of L. (This makes sense because the fibred-scattering structure 
at mf and away from sf is just the scattering structure.) The bundle S^ m \L) in the fibred- 
scattering setting is given by S^(L) = f^ /2 (L) <g> \N* mf \ m ~ N ' A ® \N* sf| m_JV/4 ®M®E, 

1 /2 

where Q b (L) denotes the 6-half-density bundle over L; a smooth nonzero section near the 
boundary is \d\dx/ x\ x ^ 2 , where A are coordinates on dL, extended into the interior of L. 
The symbol determines u G I m ' r (L) up to an element of I m+1,r (L). 



3. Geometry of the time-dependent Schrodinger operator 

3.1. Flowout from the diagonal. Let P = t 2 (D t + H). We consider this operator as an 
element of order (2, 0) of the scattering calculus for the manifold with boundary X° x X° x 
K>o- We shall use coordinates z, w, t, rj, r in the scattering cotangent bundle, as in (|2.2|l . 
The boundary symbol of P at t = is 

1 2 

P(z,w,(,V,t) = -T+-ICL- 



2 



The rescaled Hamiltonian flow is given by 



(31) 2 dC 2 dz 

i = t t = \(\l. 

The free Schrodinger propagator c„t _n / 2 e l ' z_u '' l 2t is a Legendrian distribution parametrized 
by the phase function \z — w\ 2 /2. We look for an analogue in the general case. 

Lemma 3.1. There exists a smooth Legendre submanifold L of sc T^ t=Q yX° x X° x R^o 
contained in {p — 0} that is parametrized by the phase function $(z,w) = d(z,w) 2 /2 near 
the diagonal z = w. 

Proof. Since p(z, w, d z & 7 d w $, $) = 0, the function $ parametrizes a Legendrian contained 
in {p = 0} in the region where $ is smooth, that is, within the injectivity radius. We define 
L to be the flowout along bicharacteristics from this region. L is given as a set by 

s 2 

(3.2) (z,Q =exp sg/2 (w,f}o), \r)\ = |C| = s, r = — , fj = -fj 0> s G [0,oo), (w, fj ) G S*X°. 

Since the Hamilton vector field never vanishes for s > 0, and since r — > oo as s — > oo, the 
entire flowout is smooth. □ 
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3.2. Behaviour of L near the boundary. We next analyze the behaviour of L near the 
boundary of X in the first factor. To do this, we introduce the Lie algebra of vector fields 
given by the C°°(X x X x Rj> )-span of the vector fields 

t 2 x 2 d t , tx 3 d x , tx 2 d y , tx 2 d w . 

There is a vector bundle of which these arc smooth sections forming a basis at every point. 
Let B denote the dual of this bundle; then a basis of the smooth sections of B is given 
locally near the boundary of X by 

dt dx dyi duij 
t 2 x 2 ' x 3 t ' x 2 t ' x 2 t ' 
Hence an arbitrary q 6 B may be written 3 
rn o\ k 1 Xdx udy edw 

(3 - 3) q = ^ d{ l ) + l^T + -^ + -^T- 

Then x, y, w, t, A, ji, £, K are smooth coordinates on B, up to the corner t = x = 0. We 
remark that the bundle B is a scaled cotangent bundle mixing scattering behaviour at the 
face t = with quadratic scattering behaviour at x = 0. 

There is a natural vector bundle isomorphism from SC T*(X° x X° x Rj> ) to over 1° 
given by 

(3.4) 0, w, i, C, ??, r) -» 0, uj, £, xry, x 2 r) 

in terms of the coordinates given above. Thus, the map rescales the cotangent variables by 
a power of x as we approach x = — linear rescaling for the spatial cotangent variables 
and quadratic for the temporal cotangent variable. This counteracts the growth of these 
variables as the boundary x — is approached, since we see from l|3.2|l that 77 grow linearly 
and t quadratically in s (which expected to be asymptotic to x^ 1 as s — > 00). Thus, if wc 
map L into B via the isomorphism 13. 411 . we can expect it to remain 'bounded' at x = 0. 
By abuse of notation, we denote the image of L under (|3.4|l also by L. The Legendrian L, 
now regarded as a submanifold of B, is generated by Hamiltonian flow of x 2 P starting from 
the diagonal. 

Owing to the special structure of the scattering metric g, the boundary symbol of 
x 2 P equals 

(3.5) -K+i(A 2 + ^V lMj ) 

where h lJ is a smooth function of x and y. This yields (rescaled) Hamilton equations in B 

x <>h ' ■' 

x = Xx \ = X + 2k — 1h %3 — 2~dx~^ i ^ j 

(3.6) t 2 Mi = (1 + 2A) Mi - ^, /<"/',/', 



w = 
i = t 



£ = (1 + 2AK 

k = A 2 + 2Ak + h 1 ^ jjiiUj. 

The next lemma justifies our hope that the Legendrian L stays bounded in B as we 
approach the boundary, x — > 0. 

Lemma 3.2. Along every non-trapped bicharacteristic ray inside the Legendrian L, we have 
£ — » 0, £ — » 0, A — s- -1, ju — » 0. 



^The variable /i here is not the same as the from Sections HI and 121 
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Proof. Along every nontrapped ray, we have x — > 0. Hence, there are points along the ray 
where x is arbitrarily small and where x < 0. Since x — Xx, there are points along the ray 
where x is arbitrarily small and where A < 0. 

Now we look at the equation for A. Since the ray lies inside the characteristic variety of 
x 2 P, we can substitute 2k = A 2 + f^iHj on the right hand side of Ij3.6|l . Hence 

x dhfi 

(3.7) A = A + A 2 - h l ° ■- - -^-HiHj, 

Since h lJ (x) is positive definite for small x, uniformly over dX (since dX is compact), 

x dh lj 

/'"/',/'.. - 2~dx~^ l ^ j ^ for X K e ' 

Thus, starting from a point on the ray where x < e and A < 0, we have 

x = Xx, X ^ A + A 2 if x < e. 

These equations imply that A remains negative and that x is decreasing from this point 
forwards on the ray. Moreover, since A + A 2 < for — 1 < A < 0, we see that limsup A < — 1. 
In particular, from some point on, A < —3/4. This implies that £ — > along the ray. 
Now consider the equation for fi. It implies that 

(h n mfij) = 2(1 + 2X)h lJ mfij + xX-Q—fiifJ,j. 

Since eventually A < —3/4 and x < e, positivity of /i y implies that the right hand side is 
^ — 1/2^ fiifXj, say, from some point on, so /i — > 0. Returning to the equation (|3.7|l for 
A, we see that /i — ► implies that A — > — 1, since H3.7(l with the /i terms removed has an 
attracting fixed point at A = — 1. □ 

3.3. Smoothness at the boundary. Throughout this subsection, we assume that the 
metric g is nontrapping. 

Let E(P) be the zero set of the boundary symbol of x 2 P on B P\ {t = 0}. We identify 
this set with <i sc T*(X x X°) by mapping 

k . 1. Xdx Lidy t;dw Xdx udy £dw 

x^ di l ) + ^H + ^H + ^T t0 — + — + — • 
In terms of coordinates, we map (x, y, w, K, A, fi, £) to (x, y, w, A, fi, £). This is an isomor- 
phism on S(P) since k is given by |(A,/x)| s /2 on S(P). The image of L under this 

identification is then a Lagrangian in q sc 7 1 *(x x X°) (which we shall denote L). We empha- 
size that this Lagrangian lies over the interior oflxl", not at the boundary. In terms of 
parametrizations, if $/£ is a local parametrization for L, then $ is a local parametrization 
of L. 

It turns out that L is not smooth at x — 0. To recover smoothness, we shall blow up the 
submanifold 

(3.8) S = {x = 0,X = -l,fi = 0,t = 0} 

where L meets the boundary (by Lemma f3. 211 . We shall use Melrose's notation and termi- 
nology; thus, we denote the manifold with corners obtained by blowing up the submanifold 
S by 

[ qsc T*(X xX°);S]; 

explicitly, this is the manifold with corners obtained by removing S and replacing it with its 
inward-pointing spherical normal bundle, with a natural C°° structure (see |13|L This 
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inward-pointing spherical normal bundle becomes a new boundary hypersurface of the 
blown-up space which we refer to as the 'front face' of the blowup. The lift to this blown-up 
space of the original boundary hypersurface dX x X° will be denoted sf. 

Lemma 3.3. Let (X,g) be a manifold with boundary, with a nontrapping scattering metric. 
Let S C i sc T*(X x X°) be given by i|3.8fl ■ where x is a boundary defining function such that 
\dx/x 2 \ g = 1 at dX. Then 

(3.9) X ts ee [* SC T*(X x X°); S] \ sf 

is an affine bundle over X x X° . isomorphic (but not naturally isomorphic) to the scattering 
cotangent bundle SC T*(X x X°). The Lagrangian L lifts to the space (|3.9|) to be a smooth 
manifold with boundary, such that dL is contained in the boundary of sc T*(X x X°) and 
such that L is transversal to the boundary there. 

Proof. Performing the blowup amounts to introducing the coordinates on X t s in a neigh- 
bourhood of the front face. The induced map from X t s to X x X° is a fibration in which 
each fibre is diffeomorphic to R 2 ™. To specify the affine structure we introduce the coor- 
dinates A = (A + l)/x, M — [i/x, 5 = £/x which are smooth on X t s up to x = 0. The 
affine structure is defined by deeming A, M and S to be affine coordinates on each fibre. 
This is well defined since a smooth change of variables in (x, y, w, A, /i, £) which is linear on 
each fibre of i sc T* (X x X°) induces a smooth change of variables in (x, y, w, A, M, S) which 
is affine on each fibre. (Observe that a coordinate change such as x — > x + ax 2 induces a 
translation A — > A + a + 0(x), so the coordinate changes are not in general linear on the 
fibres.) In terms of the new coordinates, a point 

X 6 x z x z 

becomes 

dx . dx dy „ dw 

(3.10) q = + A-,- + A/ ■ + H ■ — . 

X a X X X 

Mapping q to Adx/x 2 + Mdy/x + Sdx/x defines an affine bundle isomorphism between the 
blown-up space X t s defined in i|3.9[) and SC T*(X x X°). It is not, however, a canonical 
isomorphism — this is clear, since the scattering cotangent bundle has a linear structure 
which X t s lacks. We remark that if we performed the blowup at the zero section of the 
quadratic-scattering cotangent bundle, that is, at A — instead of A = —1, then the resulting 
space as in l|3.9|l would be naturally isomorphic to the scattering cotangent bundle. 

Next, we show that L has the regularity specified in the statement of the lemma. Our 
method of proof is to lift the vector field V in (231) to the space [^ SC T*(X xX°);S], and 
divide by a boundary defining function for the front face, obtaining a new vector field W. 
Clearly integral curves of V are the same as integral curves of W. We claim that W is a 
smooth vector field on [ c i sc T*(A x X°); S] which is tangent to sf but transverse to the front 
face. Consequently, every integral curve of W must meet the boundary of ^ SC T* (X x X°); S~\ 
in the interior of the front face. Moreover, the fact that W is smooth and nonvanishing at the 
front face implies that the flowout L is smooth up to the front face and meets it transversally. 
Hence to prove the lemma, it is sufficient to prove the claim above. 

Proving the claim is simply a matter of computing the vector field (|3.6() in various co- 
ordinate systems valid in different regions of the blown-up space ^ SC T*(X x X°);S\. In a 
neighbourhood of any point in the interior of the front face, we may use coordinates A, M , S 
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defined above. In terms of these coordinates, the vector field takes the form 

r dh^ 

(3.11) V=(-x + Ax 2 )d x + xW J Mjd Vi + xi-h^M.Mj + — — — M j Mj ) d\ 

1 dh? k 

+ x(AM t - — — — Mj Mk ) c?m< + zAS . d s . 

2 % 

Dividing by x, which is locally a boundary defining function for the front face in this region, 
shows that W = V/x is a smooth vector field in this region which is transverse to the 
boundary. 

Near the corner of [i sc T*(X x X°);S~\, i.e. near the boundary of the front face, we 
need several different coordinate patches. Let p = \J x 2 + (A + l) 2 + |£| 2 + \fi\ 2 , and let 
s = A + 1. When s p/4, then coordinates S\,ri = x/s%,Mi = /i/si,Si = £/si may 
be used. When |/x| ^ p/4, then valid coordinates are S2 — \p\,r2 — x/s2,A2 = (A + 
l)/s2,/t = /i/s2,S2 = £/s2 may be used, and when |£| ^ p/4, then valid coordinates are 
s 3 = lCl) r 3 = x/s^,A^ = (A + l)/s3,M3 = /i/s3, ^ = ^/s3- The union of these coordinate 
patches covers a neighbourhood of the corner of [<i sc T*(X x X°);S]. 

In the first set of coordinates, (si, r\, Mi, Hi), the vector field V takes the form (where 
the subscript 1 is dropped for notational convenience) 

(3.12) V = (s + s 2 k)d s - rskdr + sh lJ Mjd Vi 

1 dh jk 

+ s((l - sk)M t - — — — Mj Mk ) 8m i + s(l - styEidBt 

The vector field W — V/s is tangent to sf and transverse to the front face in the region of 
validity of these coordinates. The calculation for the second and third set of coordinates is 
similar, and is left to the reader. □ 

Now let G be the boundary of L: G = dL C sc T§ XxXO (X x X°). 

Lemma 3.4. The submanifold G of sc Tg XxXO (X x X°) is Legendrian. 

Proof. The contact form on sc Tg XxXO (X x X°) is given by ui(x 2 d x , •), where the symplectic 
form a; is now given by 

. . dA A dx dM Ady c?S A dw Mdx A dy Edx A dw 

(3.13) U = 2 1 ~ ^ 2 2 

X X *Xj X 

which we obtain by taking the differential of (|3.1UI) . Since V in l|3.11|l is of the form 
— xd x + xV", where V" is tangent to the boundary of SC T*(X x X°), and xuj by l|3.13|l is 
nondegenerate on vector fields tangent to the boundary, uniformly up to the boundary, we 
find that the contraction of uj and x 2 d x satisfies 

Lo(x 2 d x , •) = — xu)(V, •) + xa 

where a is a smooth one- form on SC T*(X x X°). Let W be a vector tangent to G. Then 

uj(x 2 d x , W) = -xlo{V, W) + xa(W), 

The first term vanishes since G C L, both V and W are tangent to L, and L is Lagrangian. 
The second term vanishes at x — 0. This proves that G is Legendrian. □ 
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3.4. Asymptotics of geodesic flow and the Legendrian G. The variables A and M 
have geometric meaning on the Legendrian G in terms of the asymptotic behaviour of 
geodesies. Notice that, by Ij3.2|l . the value of r on the Legendrian L is given by d(z, w) 2 /2. 
Thus, K is given by 



2 2 
K = X T = X 



d(z, w) 2 



2 

Near the boundary G of the Legendrian, we have by (|3.5|) 

K = i(A 2 + h^fiitij) = ^((-1 + xA) 2 + x 2 h tj MiMj) . 

This implies that 

(3.14) d(z,w) = - - A + 0(x). 

x 

Thus we interpret —A on G as a 'sojourn time', the time that the geodesic spends in the 
interior of the manifold before emerging into the conic region. 

Remark. It is instructive to note that under a change of boundary defining function x —> 
x = x + ax 2 which preserves the normal form l|l.ll) we find that the sojourn time changes 
as follows: letting t be a unit speed parameter along the geodesic flow, we compute 

1 ,. lx ,. 1 
limi — — = hmi — = —a + hint : 

x xx x 

hence we easily verify that this 'sojourn time,' while ill-defined as a function, does transform 
as a coordinate on the affine bundle X t s\x=0: a $ we know it must. 

To interpret the coordinate M, we notice that according to (|3.11|l we have 

//, = //'•'/'; = xh ij Mj, 

0(x 2 ), 



(3.15) ni „ which implies dyt/dx = —h^Mj + 0(x). 



Hence M measures the angle at which the geodesic strikes the boundary, or more geomet- 
rically, M is a coordinate on the pencil of 'asymptotically parallel' rays with a fixed final 
direction y. The coordinate y, of course, is just the limiting direction of the geodesic, while 
w and — S, which are constant under the flow, give the initial conditions of the geodesic. 

3.5. Fibred scattering structure. There is yet another way of viewing the Legendrian 
L, in terms of the fibred-scattering structure of X 2 described in Section |2~31 We show that 
one can map the Legendrian L C B to a fibred-scattering Legendrian of ^T*X 2 . To 
define this map, we use local coordinates (x, y, w, t). Let 

k 1 Xdx udy £dw 
(3.16) q = — d(-) + — 3- + -2- + —3-, 

x z t x a t x z t x z t 

be a point in B, using coordinates as in lEPll . We map q G L c B to 

(3-17) -"- + \-^ + S-,'--\ d ^- ,T " x -- 

Lemma 3.5. Assume that the metric g is nontrapping. Then the image of (|3.17fl is a 
smooth fibred-scattering Legendrian submanifold of^T*X 2 . 

Remark. Notice that the terms dt/t 2 x 2 and dx/x 3 t in (|3.16(l are too 'big' for the fibred 
scattering space ^T*X 2 . The addition of terms as in (|3.17|) cancels the big part of these 
terms and yields something which remains 'bounded' in <PT*Xf up to x = 0. 
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Proof. We use the blowup coordinates A, M, S as before. Thus, 
A = -1 + xA 
fi = xM 
£ = xZ 

k = i(A 2 + h^fHfj,/) = ~ - xA + x 2 (A 2 + h lj MiMj) 

Now we add \dtjx 2 t 2 + dx/x 3 t to q in l|3.16|) . and get 

. xdt + tdx , . 2 , , -,, ■ r , , . , , _ du> 

Thus there is a diffeomorphism between L and ^L, given in local coordinates by 

(s,y,«; > A,Af > a) ^ (x,y,w,A,A 2 + h ij M i M j ,M,E), 

so is smooth up to the boundary and meets it transversally. The submanifold is a 
Legendrian distribution. Indeed, if \1/ locally parametrizes L away from the boundary of B, 
then ^ — l/(2x 2 t) locally parametrizes ^L. Finally, under our identifications, the image of 
the boundary of under (3 is identical with G from Lemma 13.41 which is a Legendrian 
submanifold. This proves that *L is a fibred-scattering Legendre submanifold. □ 

3.6. Trapped rays. If the metric g has trapped rays, then we need to localize the con- 
structions of this section in the nontrapping region. We define L to be as in Lemma [3. II 
but taking only those (w, ff) £ F, rather than all (w, ff) £ S*X°. We then define L^L 
and G in terms of L as before. Then Lemmas 13.31 13.41 and 13.51 continue to hold except we 
lose smoothness of L, L and at the set corresponding to s = in (|3.2|) . Moreover, if 
we localize to a compact subset of F, then the corresponding pieces of L, L and are 
compact. 

Lemma 3.6. The map Sf defined by 

(3.18) S f (w,Z) = (y,A,M)^(y,w,A,M,-E)eG, ^ = S 

is a contact diffeomorphism S*X° D F t-> sc Tg X (X) satisfying (|1.4|) and (|1.5(l . Similarly, 
Sb{w,E) — —Sf(w,S) is a contact diffeomorphism. 

Proof. For each (w, S) € F there is a unique point (y, w, —A, — M, S) G G found by following 
the bicharacteristic emanating from (w, S) until it hits the front face of A" t g; uniqueness and 
smoothness are a consequence of the transversality of the vector field W in the proof of 
Lemma f3. 31 Invertibility of the map Sf follows from the previous discussion. In fact, given 
(y, A, M) G sc Tg X (X), y determines a pencil of geodesies with asymptotic direction y, M 
picks out a unique geodesic within the pencil and then A indicates how far along the geodesic 
one must travel to get to the initial point w. The direction S at w is the direction of that 
geodesic at w. Hence for each (y, A, M) there is a unique (w, S) with Sf(w, S) = (y, A, M). 
Comparison between (12. 1|) and l|3.10|l shows that A = — v when interpreted as a coordinate 
on SC T*X, hence 1)3.14(1 implies p.4|) . while (|3.15jl implies Ijl.Sjl . 

The fact that Sf is contact follows directly from the definition H3.18(l and the fact that 
G is Legendrian (see pl49, for a discussion in the symplectic case). □ 
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4. Propagator (nontrapping case) 

In this section, we shall construct a parametrix U{z, w, t) for the Schrodinger propagator, 
where we restrict the second variable w to some arbitrary open set Q with compact closure 
inl°. We require that U solves 

(4.1) PU{z, w, t) = t 2 {D t + H x )U(z, w, t) = £(z, w, t) G t N x N C°°(X x Q x R >0 ) for all N, 

or in other words, the error term £ is smooth and rapidly decreasing both as t — > and as 
the z variable tends to infinity. The restriction to t > is only for convenience. In addition 
U should satisfy the initial condition 

(4.2) U(-,w,t) -> 5 W in C _oc (X) as t 0, for all w G Q. 

In this section, we deal with the case in which the metric is globally nontrapping; we sketch 
the changes necessary to localize the construction in F in the following section. 

The parametrix will be a sum of four terms, ti = U\ + IA2 + U3 + U4. Correspondingly, 
the construction is divided into four steps. We first construct a U\ supported near t = 
and near the diagonal, satisfying the initial condition l|4.2|l and solving 

(4.3) PUx = Ex+nx, 

where £\ satisfies the condition for the error term in 1)4. and lZ\{z,w,t) is supported in 
i(X,g)/4 ^ d(z, w) ^ t(X,g)/2, where 1 is the injectivity radius of (X, g). In the second 
step, we solve 

(4.4) PU 2 = -11X+S2+K2, 

with zero initial conditions, where £2 satisfies the condition for the error term in 1)4. 
and with TZ2(z,w 1 t) supported in d(z,w) i{X,g)/A, t $J 1 and with x 2e. (Here, 
and henceforth, we use (x, y) as local coordinates for the z variable when it is close to the 
boundary.) In the third step, we solve 

(4.5) PU 3 = -Tl2+£ 3 +n 3 , 

with zero initial conditions, where £3 satisfies the condition for the error term in 1)4.1(1 . and 
with 7Z 3 (z, w, t) supported in x 4e. Finally, we solve 

(4.6) PU 4 = -H 3 + £4, 

with zero initial conditions, where £4 satisfies the condition for the error term in ()4.1)l . 
Clearly the sum U\ + U2 + U 3 + U4 satisfies the conditions for a parametrix. 



4.1. Step 1 — near the diagonal. We start by constructing a formal solution to 
PUi{z,w,t) = t 2 {D t + H)Ui(z,w,t) = 0, t > 0, 

near z = w, with initial condition a delta function 5 w {z). Our ansatz is 

00 

(4.7) Id = t -»/a e »(*,«)/« t j aj (z, w). 
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By a formal solution we mean that each coefficient of P vanishes (near z — w) after applying 
the operator. Applying the operator, we find 

oo 

t n l 2 PU x = e i *(*- ,0 >/ t (-* + g(V 2 $, V,$)) tja i 

j=o 

(4.8) + it ( - a(V$, V)a + |a + ^(A$)a ) 

j'=i 

This gives us a sequence of equations to be solved, one for each power oft. The first equation 
is the eikonal equation, — $ + g(V 2 $, V 2 $) = 0. This has an exact solution 

(4.9) $>= l -d{z,wf, 

which is smooth when d(z, w) is smaller than the injectivity radius of (X, g). Motivated by 
the form of the free propagator on K™, we let this be our $. 

The coefficients of P in the remainder are successively transport equations for <Zq, a%, and 
so on. Consider the transport equation for oq. Fix a w and choose normal coordinates for 
z centred at z = w. Then 

A = J^- ^ + C(z)£> z and w ) = M 2 /2 + C(|z| 3 ), 

so 

AMz,w)+n = 0(z) and 3 (V$,V) =^(z,+0(|z| 2 ))d Zj . 
Thus, the transport equation for ao has the form 

( Zi + O(\z\ 2 ))^-a = / • a , / = ^A* + J = 0(*), 
azi 2 2 

where all terms are smooth. This has a unique smooth solution satisfying 

(4.10) a Q (z, z) = 1 for all z E X°. 

By the stationary phase lemma (see for example Theorem 7.7.5), if $ and ao satisfy 
(14.91) and (|4.1(J|) and all a,j are supported within half the injectivity radius of the diagonal, 
the initial condition 1)4. 2|) is satisfied. 

The other transport equations take the form 

(zi + 0(\z\ 2 ))-^-aj +jaj = / -dj - ^Aj-Oj-i - iVaj-i, f as above. 

We inductively suppose that ao, . . . , dj—x are smooth. Then there is a unique smooth solu- 
tion aj of this equation, establishing the inductive hypothesis for dj. 

To define U\, we take our formal solution Vt±, and multiply by a smooth function x{ z i w ) 
which is equal to 1 when d(z, w) ^ l(X, o)/4, and equal to when d(z, w) ^ l(X, g)/2. Then 
we take U\ to be an asymptotic sum of the formal series so obtained. The error term we 
decompose into £ i + IZi, where 1Z\ is given by all terms containing a derivative of x, and E\ 
is the remainder. Since U\ is a formal solution, E\ is 0{t°°) as t — > 0, while 1Z\ is supported 
away from the diagonal. This completes Step 1. 
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4.2. Step 2 — at the main face. Here we solve away the error term IZi. We regard IZi 

1/2 

as a half-density by multiplication by the Riemannian half-density \dg z dg w dt\ . Then IZi 
has an asymptotic expansion of the form 

oc 

(4.11) t -n/2+i e i*{z,v,)/t i? Ti {z, w)\dg z dg w dt\ 1/2 . 

3=0 

It is therefore a Legendrian distribution, associated to the Legendrian L of Lemma 13.11 
see or section IzTTl for a brief description. 

We now need to specify how our differential operator P acts on half-densities; we shall do 
this by specifying a flat (i.e. covariant constant) half-density. One natural choice would be 

1/2 

to specify that the Riemannian half-density \dg z dg w dt\ 1 is flat. However, we shall make a 
different choice here since we want to apply the symbol calculus from 7 , where it is assumed 
that the flat half-density is a scattering density (that is, a bounded non-vanishing section 
of the scattering half-density bundle). Thus, we specify instead that 



(4.12) 



dg z dg w dt 



For convenience we define 



1/2 



is flat, i.e. V 



dg z dg w dt 



dg z dg w dt 



(4.13) 

Hence, we want to solve 
(4.14) 
or equivalently 

(4.15) (P + i(n + l)i) (K 2 \dg z dg w dt\ 1/2 



f2n+2 



1/2 



1/2 



0. 



f2n+2 



P(K 2 a) = -Kia, 



-Ui\dg z dg w dt\ 



1/2 



The additional term i(n+ l)t is a subprincipal term, since it vanishes to first order at t = 0. 

Writing 72. i in terms of a and using l|2.14|l we see that 1Z\ € I 7 / 4 (L). The Legendrian L 
is characteristic for P (that is, the symbol p of P vanishes on L, by construction of L), so 
we look for U2 £ 7 3//4 (L; ^fi 1 / 2 ) whose symbol satisfies the transport equation Ij2.7|l along 
L. In our situation, sc H p is given by l|3.1|) . while dp/dr = 1, and the subprincipal symbol 
Psub = * according to formula (2.9) of |7]. Therefore, by l|3.1|l . we can solve away the error 
term IZi G I 7 ^ 4 (L) with a Legendrian distribution u of order 3/4, by solving the ODE 



(4.16) 



l (u ) - a 7 / 4 (fti) 

with 'initial condition' that the symbol of Uq vanishes in a neighbourhood of {z — w, £ = 
77 = 0}. All bicharacteristics originate here and tend to the boundary {x = 0} by the non- 
trapping assumption, so there is a unique smooth solution with this property. It follows 
that in the region x ^ e, we have 

Piio-^i El 7/4+1 (L) 
Inductively, assuming that Uk £ I 3 ^ 4 (L) solves 

Puk-K^f'^+^L) 
in the region x e, we can find a Uk+i 6 I 3 ^ 4 (L) which solves 

Pu^-K^f/^iL) 
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in the region x e, by solving 

(4.17) -i (c Hp + k- |) a^+ k +\v k ) = a^+ k+1 (Pu k Tlx) 

and letting Uk+\ = Uk+Vk- The Vk can then be asymptotically summed to yield a«£ J 3 / 4 (L) 
which solves 

Pu-Ki<E t N C°°({x > e} x x M^ ) for all AT. 

We choose a cutoff function X2 G C°°(X) which is equal to 1 when x ^ 2e and when x ^ e. 
Define 

^2 = X2", 

then U2 satisfies 

PU 2 = -K^Ki+Si, 

where IZ2 is a Legendre distribution of order 7/4 supported in e ^ x ^ 2e (the 7?-2 error 
comes from derivatives hitting the cutoff function x%)- Our new error term IZ2 is now 
localized near the corner. 

4.3. Step 3 — near the corner. Recall that we expect the fundamental solution to have 
quadratic oscillations at spatial infinity, i.e. at x = 0. In fact, we expect that e~ l ^ 2tx U(t) 
is a somewhat simpler kernel than U itself, near infinity. Therefore, we look for U% in the 
form 

U 3 \dg z dg w dt\ 1/2 = e^ 2tx2 U 3 \dg z dg w dt\ 1/2 , 

where U 3 solves 

Here a is the half-density i|4.13[l . Then in terms of U 3 we have 

(2 Tic \ ~ / 2 

t 2 D t - txD x + t 2 H + i(n + l)t - — JW 3 = - e - t/2tx TZ 2 . 

Here the term i(n + l)t comes from the fact that we have chosen the scattering half- 
density a rather than the Riemannian half-density to be covariantly constant, as in Step 
2. Note that the differential operators appearing in (|4.18(l are generated by the fibred- 
scattering vector fields on X 2 appearing in l|2.11|l . This is not too surprising, since the 
analytic operation of multiplying U 3 by e~ l l 2tx2 corresponds to the geometric operation of 
(|3.17l) which maps the Lagrangian L to a smooth Legendrian in ^T*X 2 . Since is 
well-behaved in ^T*X 2 , in the sense of being smooth up to the boundary, we can expect 
the corresponding operator (|4.18|) to be well-behaved in terms of the fibred-scatteringfibred- 
scattering structure on !, 'T*X 2 . 

Lemma 4.1. The kernel e ~ l / 2tx ' 2 7?. 2 is an element of I 7 / 4 '°°(^L), as defined in 0. 

Proof. If ^ parametrizes the Legendrian L, then W — l/(2tx 2 ) parametrizes according 
to Lemma 13.51 Hence e~ l ^ 2tx 1Z 2 is Legendrian with respect to ^L, of order 7/4 at mf 
according to the calculation in Step 2. Since H.2 vanishes in a neighbourhood of x = 0, it is 
order 00 at sf. □ 

Thus, it makes sense to seek the solution U 3 to l|4.18[l in the space pi/4,* (4 L). Here 
the value of *, specifying the rate of decay of the symbol as x — > 0, is determined by the 
transport equation which is regular singular at x = 0. 
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To do this, we must analyze the Hamilton vector field of our operator on near the 
boundary at x = 0. This is given (modulo some identifications) by the vector field V in 
(|3.11|1 , which is xd x + xW, where W is tangent to the boundary. To determine the sub- 
principal symbol, observe that the adjoint of t(xD x — in/2) with respect to the Riemannian 
density is t(xD x — in/2) + 0(tx) near the corner, hence the Weyl symbol of t(xD x — in/2) is 
v + 0(tx) where v is the variable dual to x. Hence the subprincipal symbol of t(xD x — in/2) 
at mf vanishes at the boundary of mf (that is, at x = 0). Also, the subprincipal symbol 
of t 2 H vanishes identically on the main face. So the boundary subprincipal symbol of the 
operator in (|4.18l) is equal to i + 0(x). Thus, we look for a u 6 7 3 / 4 >*(*L) satisfying 

(4.19) ( - iC- xdx+a: w ~ i{\ + \ ~ ^ 1 ) + * + xq>)* 3 / 4 (u ) = a 7 / 4 (e-^ tx2 TZ 2 ) 
This gives an equation of the form 

(Cza.+*w + | + xq')a 3 / 4 (u ) = -iaV\e- i ' 2b *'R 2 ), 

so by (|2.6|) . cr 3 / 4 (ito) is of the form aox~ n / 2 \x~ 1 dxd\\ 1 ^ 2 \dt\ 3 ^ 4 ( 2n+1 )/ 4 j w here ao is smooth, 
and A are coordinates on the boundary of ^L, extended into the interior. Here we have 
explicitly included the power of \dt\ in the formula since we now need to change boundary 
defining function from t to p — tx, in order to apply formulae from [7|. Doing this, we find 

that <7 3 / 4 (iio) is of the form aQX~ 1 / 2 \x~ 1 dxd\^ 2 \dp^' i ^^ 2n+1 ^ A , where ao is smooth. By 
the (corrected) symbol calculus for fibred Legendrians given by Proposition 3.4 of [7] , this 
implies that u e J 3 / 4 - 1 / 4 ^), and that 

(4.20) (t 2 D t - txD x + t 2 H + i(n + l)t - ^)u + e^^Tli G f^+^^L). 
Inductively, we look for ui e / 3 /4,i/4(* jL ) such that 

(4.21) (t 2 B t - txD x + t 2 H + i(n + l)i - ^)«, + e-'^Tl^ G /W+M/y ^ 

In fact, we will show more — we will show that there is ui as above whose (full) symbol is 
of the form 

(4.22) a{p,x)x- 1/2 \x- 1 dxd\\ 1/2 \d P \ 3/4 ~ (2n+1)/ \ with a smooth. 

The important point here is that a is a smooth function of p and x, not just a smooth 
function of t and x. To show this, we use the boundary defining function p rather than t 
for the main face, even though it is degenerate at the corner t = x = 0. In terms of the 
coordinates p, x, y, w near the corner, partial derivatives transform as 





tD t \ 


x,y 


= PD P \ 




(4.23) 


XDy 


\x,t 


= XDy 


1 p,x 




xD x 


\y,t 


= xD x 


+ pD p \ 

\y,p r H lx,y 



We also write V = x 2 V where V G C°°(X) by assumption. Hence the operator in (|4.21|) 
takes the form 

o 2 / ~ \ n + 2 

(4.24) P = - P D X + V— [{xD x + pD p ) 2 + i(n - 2)(xD x + P D P ) + A, l(x) + V) + W^r- 

4 The power of p in the final nonzero term of the exact sequence of Proposition 3.4 of |7| is incorrect. It 
should be r — m, not m — r. 
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Any Legendrian distribution u in / 3 / 4 ' 1 / 4 (*L) can be written, modulo x°°t 00 C°°(X x K^o), 
as a finite sum of oscillatory integrals of the form (written as a Taylor series in p) 

p oo 

(4.25) u= e i ^ x ' y > w > v V p p n / 2+1 - k / 2 Y'p j a j (x,y,w,v)dv 

where ip parametrizes ^L, K C K fc is bounded, and a, £ x - ™/ 2-1- ^ 00 . We claim that in 
fact, a, 6 a; - "/ 2-1 ^ 30 , or in other words, that the full symbol of u is smooth in x and p. 
The equation to be solved is 

p OO 

(4.26) Pu = -e' l/2tx2 TZ 2 = / e^ (x ' y ' w ' v)/p p n/2+2 - k/2 J2 f^b^x^y^w^v) dv 

where the bj are supported in x ^ e. Substituting (|4.24() for P, we find that the left hand 

side is given by 

(4.27) 



e 

K 



i^,y, W , v)/Pp n /2+ l- k/ 2\ ( _ ^ + 1 ( ^ _ + + p (_ Dx + 



{xipx ~ ip){xD x + pD p ) + (V y ip,D y ) h ( x) + v-- h/J + p 2 Q \ y^ j p i a i {x,y,w,v)dv, 

X J j=o 

where Q is a second order operator generated over C°° by ie-D^, pZ? p , and D y . Since ip 
parametrizes a Legendrian which is characteristic for P, the expression 

(4.28) -^x + ~(^*-V0 2 + ~MIm*) 

vanishes when d v tp = 0. As ip satisfies the nondegeneracy condition (|2.13|) . the function 
(|4.28|) can be expressed 

E" dip . , , 
a i~R — with ai smooth. 
OVl 

i=i 

We then write ai(d Vi ip)e l ^/ p = aipD Vi e^' ' p and integrate by parts in v. These terms then 
become 0{p) terms. Let us write W = (xip x — ip)xD x + (V y ip, Dy)^/^ + Y2i a iD Vl - The 
equation becomes 

oo ( 2 ^ oo 

(4.29) Y,P j \(-D x + W + i^— + /,-) + P Q, \a 3 =J2 A • 

3=0 { ) j=0 

Here, / and fj are smooth and Qj is a second order operator generated by xD Xl D y and D v 
with smooth coefficients. Thus, we need to solve 

n + 2 

{-xD x + xW + i— h xf )a = xb Q , 

( 4 ' 3 °) n + 2 

{—xD x + xW + i — h xfj)aj = xbj + xQjdj-i. 

Here bj is supported in x £ [e, 2e] and our initial condition for a,j is that it vanish for x > 2e. 
These are regular singular ODEs for the a,j. The first equation certainly has a solution 
which is in x~ n / 2 ~ 1 C oc . Assuming that a 3 is in x~ n t 2 ~ l C°° , it follows from l|4.30() that 
a J+ i S x~ n l 2 ~ 1 G oa , since the right hand side of l|4.30() is in x~ n l 2 C° a and the indicial root 
of the regular singular operator on the left hand side is — n/2 — 1. Hence the claim follows 
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by induction. (Notice that this induction would not work if we did not make the stronger 
assumption 14.22fl on the symbol of u.) 

Thus, we let Uz be an asymptotic sum of the formal series (|4.25(l constructed above, 
supported say in t ^ 1. This is certainly in / 3 / 4 ' 1 / 4 ( ( ^L) and solves the equation 

t 2 D t - txD x + t 2 H + i(n + l)t - —)U Z = - e - l/2tx TZ 2 + £ 3 + R<3, 

where TZ 3 £ I 00 ' 1 / 4 This completes Step 3. 

4.4. Step 4 — near spatial infinity. Now we have an error term 71$ S I 00,1 / 4 ^!/), and 
we seek a solution to 

(4.32) (t 2 D t - txD x + t 2 H + it^^JUi = -TZ 3 € x^t^C^iX 2 ). 

We seek a solution U4 in the space / 0O ' 1 / 4 (^L). 

We can write —IZ3 as a finite sum of terms of the form 
« 00 

(4.33) / e l ^ x > y > w > v)/p p~ k/2 ^r P >b J {t,y,w,v)dv 
Jk 3=0 

where we have expanded the symbol as a Taylor series in p. Here each bj is 0(t°°) at t = 0. 
In terms of the coordinates (t, p, y), 

tD t \ =tD t \ -xD x \ 4 , 

z \p,y z \x,y ■ L \t,y' 

so the operator in l|4.32|l takes the form 

(4.34) t 2 D t \ p + it^p + ( P 2 D P ) 2 + W{pD Vi ){pD Vj ) + 

where Q is a second order differential operator generated over C°° by the vector fields p 2 D p 
and pD y . Hence, we wish to solve 
(4-35) " 

(t 2 D t \ p + it^ + { P 2 D p f + Wp 2 D Vi D yj +&)[ e^'wU" £ p>^a 3 (t, y, w, v) dv 

K J=° 



This gives us equations 

. n + 2 



( 2 n + 1 \ 
\t D t + i<— 1" c \Vi w,v)ja = b a , 



(4 - 36) / n + 2 \ 

[t 2 D t + it— h c{y,w,v)^a,j = bj + Rj{a , . . .,aj-i), 

where c(y,w,v) = ijj(0,y,w,v) 2 + h ls ipy i (0, y, w, v)ip y , (0, y, w, v) is independent of t. The 
operator Rj is such that P +1 Rj is a second order differential operator in y with smooth 
coefficients. We show inductively that there is a solution with each a,j £ t°°C°°. Indeed, 
the equation for ao is explicitly solved by 

ao(t) = f(«+ 2 )/ 2 e ic /* f ie ~^ s bo(s) S ~ n/2+3 ds 
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which is in t^C 00 since 60 G t°°C°° . Inductively assuming that do, • • • , dj— 1 are in t^C 00 , 
then it follows from the formula above, with 60 replaced by &j + Rj(ao, . . . , aj-i) that 
Oj € i^C 00 . Let be an asymptotic sum of the series on the left hand side of l|4.35l) . 
Then solves <|4.32[1 up to an x^t^C 00 error term. This completes the construction of 
the parametrix in the nontrapping case. 



4.5. Exact solution. We have now shown that 

(4.37) {D t + H)U{x,y,w,t) = e{x,y,w,t) G C 00 ^ x X x G) = C°°(M >0 ; C°°(X x £)). 
The exact propagator U(t) is given in terms of the parametrix lA(i) by Duhamel's formula 

(4.38) U(t)=K{t) + i [ U(s)(D t + H)U(t-s)ds. 

Jo 

It follows from a commutator argument due to Craig, J2J Theoreme 14, that, letting 

k 

s=0 

we have for all k > 

(4.39) U(t):H k ^L£ c (R t ;H k ). 
Note that f] k Hk — C°°, so in particular 

(4.40) U(t) : C°°(X) -> C°°(X) for all t 0. 
Lemma 4.2. Let e G C°°{R >Q ; C°°{X x G)). Then 

Ke(t)(z,w)= f ds f U{s){z,z')e{t-s){z',w)dg{z')&C oo {^ -C oo {Xxg)). 
Jo Jx 

Proof. Let us fix w G Q; then we may regard e(z,w,t) as an element of £"C°°(R^o; "Hk) for 
any ra, k G N. First we establish that #e G W°(l ?0 ; H k ). We compute for i G [0, T] 

ll^e(t)|| w> < / ||E/( S )e(t- S )|| w 







s^C T / ||e(i-a 
Jo 

= 0{t n ) for all n G N 



This shows that Ke G t n L°°(t; Hk) for every A: and n. Derivatives of t and w can now 
be estimated similarly and uniform estimates then follow from the relative compactness of 

Q. □ 

We conclude from pg| and Lemma IO that U(t) -U(t) G C°° (R >0 ; C°° (X x (?))■ This 
proves that the kernel of U(t) is such that e" i/2a;2t C/(i) G / 3/4 ' 1 / 4 ( <#, L), which completes the 
proof of Theorem ll.il 
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5. The trapping case 

We now sketch briefly the changes to be made in the parametrix construction in case 
there exist trapped geodesies. In this case, given a properly supported Z £ ^ a (G) with 
WF'Z C F, the non-forward-trapped part of the phase space, we wish to construct a partial 
parametrix Uz such that 

(D t + H)U z (z,w,t) = e(z,w,t) G C°° (R >0 ; C°° (X x Q)). 

with initial condition 

U z {t)f -> Zf in C-°°{X) as t -> 0. 

It then follows by Duhamel's principle that U z - U{t)Z G C' X {R >0 ;C°°{X x Q)). 

The only modification necessary in the construction comes at the end of Step 1. We 
begin with a result about the composition of pseudodifferential operators in the boundary 
with the simple Legendrian distributions appearing in Step 1. 

Lemma 5.1. Let k(z,w,t) — t s a(z,w,t)e i ^^ z ' w ^' t such that a and <I> are smooth, a is 
compactly supported and d w $ / on the support of a. Let B G (X) be a properly 
supported, classical pseudodifferential operator. Then 

v = k(z,w,t)B(w,w') dg(w) 
Jx 

has the form 

v = t s a(z, w, t)e i<Sf <- z ' w V t , a smooth, 

with 

WF S » C WF sc (k) n {(z, w, C, rj, t) I (w, 77) G WF'(B)}. 

Proof. The result follows directly from the lemma of stationary phase. Alternatively, we 
may regard z and t as smooth parameters and regard v as B T (k), the transpose of B acting 
on the distribution k. Recall that when B is a zeroth order classical pseudodifferential 
operator and c,(f> G C°° with c compactly supported and d(p ^ on supp c, 

B{c{w)e^ [w)lt ) =c(w,t)e^ [w)jt 

where c G C°°(R^ x X) and c(w,t) = 0(t°°) for all w such that (w,d w (f>) (£ WF'B. (See 
Chapter 8, section 7.) The lemma now follows directly. □ 

To construct Uz, we construct U\ as in Step 1 of the previous section and consider IA\ o Z . 
The error term o Z = (D t + H)Ui o Z is then a Legendre distribution associated to L, by 
the lemma just proved (since d w & ^ on the support of the symbol of IZi). Moreover the 
lemma shows that the scattering wavefront set of (D t + H)U\ o Z lies on the portion of L 
emanating from (w, rf) G WF'(Z), which is contained in the non-forward trapped set. Now 
we proceed with Steps 2, 3 and 4 with TZ\ o Z substituted for TZ\. As no trapped rays lie 
in the microsupport of 1Z\ o Z however, we need only solve the transport equation in Step 
2 along the remaining, non-trapped, rays. The rest of the construction proceeds as before. 

Thus, we have shown that Uz is the sum of two terms, one of which has the form U\o Z 
where U\ is a Legendrian distribution as in subsection l4.ll and the other is a fibred-scattering 
Legendrian distribution of order (|, j). 

Remark. It is not true that U\ o Z is a Legendrian distribution associated to L, since 
d z & = at z — w. In fact, U\ o Z is a Legendrian conic pair associated to (A, L), where 
A is the diagonal Legendrian {(z, £, z, — £, k = 0)}. Since we do not need this fact, we omit 
the proof. 
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6. EGOROV THEOREM 

Let W(t) — e~ l l 2tx2 ~U{t). Suppose A is a properly supported pseudodifferential operator 
on X°, microsupported inOxO with O C B. We will show that B = W(t)AW(t)* is 
a scattering pseudodifferential operator of the same order as A, and that the symbol and 
microsupport of B are given by the pullback of the symbol and microsupport of A by the 
sojourn relation, scaled by the factor t . Indeed we prove a more general result. To state it, 
suppose that S is a contact diffeomorphism from O C S*X° — > sc Tg X X. Then, by definition, 
S*x, the pullback of the contact form on sc Tq X X, is equal to a multiple ffj-dw of the contact 
form on S*X. Associated to S is a Legendrian submanifold G of sc Tg XxXO X x X°, given 
by 

(6.1) G = {(y, w, v, fi, -frj) | (w, rj) G O, (y, v, fi) = S(w, i])}. 
Indeed this is the correspondence of Lemma |3~B1 

Proposition 6.1. Let S be a contact diffeomorphism S : O C S*X° — > sc Tg X X and let G 
be the Legendre submanifold given by (EU. Suppose that W G I°{X x X°,G; sc f^ x Qti), 
and suppose that A £ ^ m (X°; fi 5 ) is properly supported and microsupported in O. Then 

(6.2) WAW* G 8c *-°°' m (X; sc n5) mith symbol \<r(W)(q G )\ 2 ■ <j(A)(S- 1 (q)) 

where qo is the point on G corresponding to q. Conversely, if A G sc vI/ _00 ' m (X; sc f2^ ^ /j as 
compact microsupport contained in the range of S, then modulo a kernel in C°°(X x X), 
W* AW G ty m (X°) is a properly supported pseudodifferential operator with symbol 

(6.3) o-(W*AW)(() = o-(A)(S(0) ■ k(^)(C G )| 2 . 

The proof proceeds as follows: we reduce to the case X = MP, the radial compactification 
of Euclidean space, by localization; applying the Fourier transform and results of Melrose- 
Zworski I14 | , we then deduce the result from the standard Egorov theorem. 

To begin, we need a lemma, in effect a version of Proposition 10 of |14j with parameters, 
which tells us what happens when we pull back a Legendrian distribution on R" x R n by the 
Fourier transform T . The Fourier transform invariantly maps half-densities on a vector space 
to half-densities on the dual space. We recall from that it can be interpreted as a Fourier 
integral operator associated to the Legendre diffeomorphism L : sc T* a jp-R" — > S'*((IR n )*) 
given by L(z, C) = {(, —z), where z G 5 n_1 , the radial compactification of R n , and Q G (R™)*, 
the dual space. 

Lemma 6.2. Let S, G and W be as in Proposition \^\ for X — R™, and let T be the 
Fourier transform. Then 

To W G L°((R n )* x m. n ,G L ;&) +<S((R")* x R"), 

where Gl is the Lagrangian associated to the graph of the contact transformation L o S, L 
is the Legendre diffeomorphism, and S denotes Schwartz space. 

Proof of Lemma \b.<H Without loss of generality we may assume that W can be written in 
terms of a single parametrization. Then W can be written 

fe /*_£./*> \ l 1 i k 1 

(27r) _ 2 / e l(p{z ' w ' v '' z 'a(z, —,w 1 v)\z\2 dv \dz dw\z , a smooth, compactly supported, 
Jr» \ z \ 

modulo S(R n x R"), so F o W is given by 

(27r) -„/2-fc/2 f I e-<«-C+<*(*.«'.*)Wo(5,|«|- 1 ,«; ) «)|2| fc / 2 d»tte|<iC|*|dt«|* 
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modulo <S((R ra )* x K n ). Imitating the proof of Proposition 10 of |TIj, we set ^ — v \ z \ so t na t 
To W is given, modulo 5((R n )* x M"), by 

(6.4, P .)^//^^«l.^, W -.^ S /|,|)|,|-V.* (bM Cl»|*.|». 

The phase is now homogeneous of degree one in (z, w) and parametrizes the Lagrangian Gl 
while the amplitude a(z, \z\~ 1 , w, v/\z\)\z\~~ fe / 2 has order —k/2 with k + n phase variables 
z, u. Hence ljfi.41) is a Lagrangian distribution of Lagrangian order 0. □ 

Proof of Provosition \6.1\ . By a microlocal partition of unity, we may assume that WF'A C 
O, a small subset of S*X°. Then S(0) lies over a small closed set in dX. Without loss of 
generality, we may assume that the support of the kernel W(z, w) in the z variable lies within 
a small neighborhood W, covered by a single coordinate chart (x, y), of this closed set, since 
introducing a spatial cutoff introduces only residual terms in WAW*. By identifying W 
with a neighbourhood of a point of the boundary of R ra and identifying a neighbourhood of 
7r(0) with a neighbourhood of the origin in M™, we may assume that that A G ^(K™) is 
properly supported and that W is as in Lemma \d. 21 
Now we write 

WAW* = T*{TW)A{TW)*T, 

and apply the standard Egorov theorem to the middle three factors to conclude that this 
takes the form 

F*AT, 

where A € * m (K™) + S(R n x R n ) is, up to a Schwartz kernel, a pseudodifferential operator 
with principal symbol 

a(q) = a(A)(L-\S-\q))) ■ a{{TW){TW)*){q) = a{A){L^ 1 {S^ 1 (q))) ■ \a{W){q G )\ 2 . 

Proposition 8 of ^1] now shows that 

T*AT E sc *-°°' m (l^) 

with principal symbol 

a{A){S-\ q ))-\a{W){ qG )\ 2 . 
This proves (|6.2f) . The proof of (|6.3|l is similar, proceeding from the expression 

W*AW = W*T*TAT*TW. 

□ 

7. Wavefront set bound 

In this section we prove Theorem II .21 First we prove a preliminary result on the propa- 
gator at a fixed time t ^ 0. We define the contact transformations St for t ^ as follows: 
for t > and q £ S*X°, let S t (q) with domain F be given by S t (q) = t~ 1 Sf(q) where Sf 
is the contact transformation defined in Lemma \'S . 61 and the scaling by t^ 1 acts in the fibre 
variables. For t < 0, let St with domain B be given by St(q) = \t\~ 1 Sb(q). Let G t , t ^ 0, 
be the Legendrian (|6.1fl determined by S t - 

Lemma 7.1. Let Z G ^°(X°) be properly supported, and micro supported inside F. For 
fixed t > 0, the operator W(t)Z is in I°(X x X°; G t ). Similarly, if Z' G is properly 

supported, and micro supported inside B, then for fixed t < 0, the operator W(t)Z' is in 
I°(X x X°;G t ). 
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Proof. The construction of the parametrix shows that W(t)Z is given by a finite sum of 
oscillatory integrals 

e iH0,y, w ,v)/ x t x -k/2 a ^ y) Wj Vj t) d , dx dy dw , 1/2 



xn +l I ' 
A'gR fe 

where ip and a are constructed in sections 14.31 and l4~4l (Here, because t > 0, we may replace 
y, if, i>) by ip(0,y,w,v) in the phase.) According to the remark below equation l|2.15|l . 
this is a Legendre distribution associated to the Legendre submanifold Gt, of order zero. 
The proof for t < is similar. □ 

The proof of Theorem ll.2l is now straightforward. 

Proof of Theorem M . 6 <A We only prove the theorem for t > 0, since the argument is similar 
for t < 0. Let W(t) = e~ l/2tx2 U(t) as above. Then since both U(t) and the multiplication 
operator e ~ l / 2tx are unitary and map C°°(X) to itself, the same is true of W(t). For 
simplicity we shall first assume that the metric is nontrapping. 

Suppose that q ^ WF(u(-,i)) for some t > 0. Then there is a properly supported 
A e ^"(X) which is elliptic at q and such that Au(-,t) £ C°°(X). 

Since / = U{-t)u, we have / = e - l / 2tx 2 \Y(-t)u, or u = W{-t)*e^ 2tx2 f . Hence 

A?i = AW{-t)*e l / 2tx2 f 6 C°°(X). 

Applying we see that also W(-t)AW{-t)*e l / 2tx2 f e C°°(X). By Proposition O 

W(— t)AW(— t)* — A is a scattering pseudodifferential operator of order zero, which is 
elliptic at S*_ t (g). We have just shown that A(e l / 2tx f) £ C°°(X), so by definition of the 
scattering wavefront set, we see that S-t(q) ^ WF sc (e l ^ 2tx f). Conversely, assume that 
S- t (q) WF sc (e'/ 2tI /). Then there exists a scattering pseudodifferential operator A with 
compact microsupport which is elliptic at S~t(q) and such that A(e l l 2tx f) £ C°°(X). By 
PropositionKHl we may write A = W(-t)AW(-t)* , with A = A' + A" S +C°°(X x 

X), where A' is properly supported and elliptic at q and A" has a Schwartz kernel. Hence 
PF(-t)AVF(~i)*e i / 2t:,;2 / = W(~t)Au(-,t) is in C°°(X), which implies (using (|PUI) ) that 
Au(-,t) G C°°(X). Hence g i WF(u(-,t)). 

Now we do not assume that the metric is nontrapping, but assume instead that the point 
q is not backward-trapped. Then since the set B of non-backward-trapped points is open, 
one can choose A as above so that WF'A C B. In fact, one can find in addition a properly 
supported pseudodifferential operator Z such that Z is equal to the identity microlocally on 
WF'(A) and such that WF' Z C B. Then ZAZu(-,t) e C°°{X). Repeating the argument 
above, we see that 

W(-t)ZAZW{-t)*e^ 2tx2 f S C°°(X). 
In the trapping case, W(—t)Z is a Legendre distribution associated to SL t , so by Proposi- 
tion l6.ll W(—t)ZAZW(—t)* is a scattering pseudodifferential operator which is elliptic at 
S-t(q), and the rest of the argument goes as before. □ 
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